Evidence of strong higher twist effects in diffractive DIS at HERA at moderate Q 2 
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We study a twist decomposition of diffractive structure functions in the diffractive deep inelastic 
scattering (DD1S) at HERA. At low Q 2 and at large energy the data exhibit a strong excess, up to 
about 100%, above the twist-two NLO DGLAP description. The excess in consistent with higher 
twist effects. It is found, that complementing the DGLAP fit by twist 4 and 6 components of the 
GBW saturation model leads to a good description of data at low Q 2 . We conclude that the DDIS 
at HERA provides the first, strong evidence of higher twist effects in DIS. 
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1. Introduction and main results. The fundamen- 
tal description of hadron scattering processes is based on 
the Wilson Operator Product Expansion within Quan- 
tum Chromodynamics (QCD). In this procedure the scat- 
tering amplitudes are expanded into a series of contri- 
butions labelled by twist r. The twist determines the 
exponent of leading order scaling behaviour of those con- 
tributions in the inverse powers of the characteristic large 
scale. In Deeply Inelastic Scattering (DIS) the large scale 
is set by the negative momentum transfer, Q 2 , from the 
electron, e, to the proton, p, corresponding to the virtu- 
ally of an exchanged photon, 7* . The twist expansion of 
the DIS cross-section takes the form a(Q 2 ) = J2t °v(Q 2 ), 
with <r T (Q 2 ) ~ 1/Q r , with r = 2,4, . ... In the leading 
contribution at large Q 2 , the proton structure enters at 
t = 2 in the form of universal parton distribution func- 
tions (pdfs), whose QCD dependence on the hard scale is 
described by DGLAP evolution equation. Although the 
twist-2 description of scattering processes is successful, 
it has important limitations, that come from neglecting 
the higher twist (HT) terms. 

In a complete and accurate description of hadron scat- 
tering the HT contribution should be taken into account. 
First, they correspond to fundamental hadronic matrix 
elements that provide additional information on hadron 
structure, beyond parton densities. Despite the power 
suppression, the HT effects become relevant below some 
scale, that depends on the process and the required pre- 
cision. In particular, a well known problem in precision 
determination of pdfs is the optimal choice of the lowest 
scale of fitted data sets. Finally, the HT contributions 
are related to multiple scattering processes, vital for good 
description of event properties at hadronic colliders. 

Thus far experimental information about HT effects in 
inclusive hadronic processes is rather poor. A number of 
phenomenological data analyses were performed for DIS 
at HERA, showing very weak HT effects for the best 
measured structure function F^ down to Q 2 — 1 GeV 2 Q] . 
For the longitudinal structure function, Fl, however, the 
HT effect are expected to become relevant below Q 2 — 
10 GeV 2 . Unfortunately, present F L data accuracy at 
low Q 2 is not sufficient for determination of HT effects. 

The diffractive DIS (DDIS) is a semi-inclusive DIS pro- 
cess in which the proton scatters elastically. DDIS events 



are quite important, as they make up to about 10% of 
DIS events at HERA. The standard QCD description 
of DDIS is based on the leading twist approximation in 
which the diffractive proton structure functions, F% and 
FjP , are expressed in terms of diffractive parton densi- 
ties. This approach is justified by Collins factorisation 
theorem for DDIS [2J. However, HT effects are not well 
understood which may lead to significant uncertainties 
and contamination in determination of diffractive pdfs. 
In particular, it is not clear how to choose the lower limit 
of the scale in DGLAP fits of diffractive pdfs. 

In this paper we perform the first dedicated study of 
HT effects in diffractive structure functions. We study a 
spectacular break down of the standard twist-2 DGLAP 
description of ZEUS data on the DDIS cross-section be- 
low Q 2 = 5 GeV 2 [31. In terms of X 2 , the DGLAP fits 
fail badly below this scale (the relative deviation reaches 
100%) and cannot be extended to the low Q 2 region. In 
the region of low Q 2 the data deviations from extrap- 
olated DGLAP predictions grow rapidly with decreas- 
ing Q 2 and with increasing photon-proton collision en- 
ergy W [J] . We find that relative deviations of the data 
from DGLAP predictions exhibit very steep, power like, 
1/Q 2 -f- l/Q 4 dependence, also the energy dependence is 
unusually steep. Such behavior of the cross-section is 
characteristic and natural for HT effects and, by far, too 
rapid for twist-2 DGLAP evolution. 

In order to estimate the HT contributions we used 
GBW saturation model @|. The model is deeply rooted 
in QCD and offers universal description of variety of ob- 
servables in DIS, DDIS and vector meson production, see 
e.g. (4j-[6]. Moreover, the structure of the GBW model 
allows for a consistent twist expansion [I]. We decom- 
posed the saturation model predictions for DDIS struc- 
ture functions into twist components and compared the 
results to the DGLAP fits and to the data. We found: (i) 
full consistence of twist-2 component of the GBW predic- 
tion with the standard DGLAP prediction — both unable 
to describe the low Q 2 DDIS data; (ii) improved, but still 
unsatisfactory description of low-Q 2 DDIS data with the 
complete (all-twists) GBW model; (iii) a good descrip- 
tion of the DDIS data with the GBW model truncated 
to include only its twist-2, 4 and 6 components. This 
success of the twist-truncated GBW model confirms that 



the data deviation from twist-2 DGLAP description are 
consistent with HT effects. We use the twist truncation 
of GBW to analyse and parameterise the deviation from 
DGLAP description of the DDIS. In fact, such a trunca- 
tion may be motivated within QCD, see the last section. 

To summarise the main results of the paper: the avail- 
able data on low Q 2 diffractive DIS show strong break- 
down of twist-2 DGLAP approximation and the devia- 
tions are consistent with sizeable twist-4 and twist-6 ef- 
fects. Thus, the DDIS data at low Q 2 provide the first 
and strong evidence for HT effects in DIS in the pertur- 
bative domain. This opens possibility for further exper- 
imental and theory investigations of the HT effects and 
for significant improvement of the diffractive pdfs. 

More details of all the calculations and results will be 
given in a coming, more extended paper. 

2. The data — the breakdown of twist-2 descrip- 
tion. The DDIS is an ep scattering process e(p e ) p(P) — > 
e(p' e )p(P')X(Px), with four-momenta p e , p' e (P, P') for 
the scattering electron (proton). The final hadronic state, 
A, with four- momentum Px, is separated in rapidity 
from the proton, that scatters elastically. The process 
is mediated by a virtual photon exchange, 7* (q), with 
q = p' e — p e . The DDIS differential cross-section is ex- 
pressed by invariants: y = (p e q)/{p e P), Q 2 = —q 2 , 
i = {Q 2 + M 2 X )/(W 2 + Q 2 ), p = Q 2 /(Q 2 + M 2 X ), and 
t = (P 1 — P) 2 , where W 2 = (p + q) 2 is the invariant mass 
squared in photon-proton scattering, and M x = P x 
is the invariant mass of the hadronic state X. The t- 
integrated ep cross-section reads: 
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where the "reduced cross-section" may be expressed 
through diffractive structure functions, F 2 and F L : 
<r?V>V3,Q*,S) = if (3) - y 2 /[l + (1 - y) 2 ]F^\ with 
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and the structure functions 



T, L may be, respectively, expressed through transver- 
saly and longitudinally polarised 7* -proton cross sec- 
tions, F^ } = (Q 4 /47r 2 a cm /3£) dal'£/dM x . 

In the recent analysis [3] the ZEUS diffractive data, 
currently the most accurate in the low Q 2 region, were fit- 
ted within NLO DGLAP approximation. The data cover 
the region of 2 < Q 2 < 305 GeV 2 but a satisfactory 
DGLAP description was found only for Q 2 > Q 2 nin = 
5 GeV with x 2 /d.o.f. rapidly growing with decreasing 
<3 2 nin . Following Ref. [3] we have calculated x 2 /d.o.f. for 
subsets of ZEUS LRG data with Q 2 > Q 2 ^, this time 
restricted to j3 > 0.035 in order to cut off contributions 
of highly resolved projectiles. We have found basically 
the same behaviour — x 2 /d.o.f. ~ 1 for Q 2 nin > 5 GeV 2 , 
and x 2 /d.o.f. reaching ~ 3 for the full Q 2 range, see the 
continuous curve in Fig. 1. 

In order to analyse the origin of the DGLAP fit dete- 
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FIG. 1. The x7<ho-i- for NLO DGLAP and NLO DGLAP 
+ HT fits to ZEUS LRG data [3] with Q 2 < Q 2 ^. 



rioration we focus on the low Q 2 region. The ZEUS fits 
[3] were performed in the 'safe' region, Q 2 > 5 GeV 2 , 
and the predictions were extrapolated to lower Q 2 . The 
deviations are found to grow rapidly with decreasing £ 
and Q 2 (see Figs. 3,5 of Ref. [3] and Fig. 2). The rel- 
ative deviation is largest, about 100%, at the minimal 
£ ~ 4 • 10~ 4 and Q 2 = 2.5 GeV 2 . 

We performed a more detailed, numerical analysis, 
assuming contributions of HTs (the details arc presented 
in the next part of the paper and Fig. 2) and we found 
that the deviation is consistent with a sum of twist-4 
and twist-6 contributions. This, in turn, implies a 
very strong dependence of the data deviation from 
twist-2 structure functions, that scales as l/Q 2 -j- 1/Q 4 - 
However, at twist 2, a possible Q 2 modification of the 
leading Q 2 dependence, ~ const (Q 2 ), due to the DGLAP 
evolution is moderate ~ Q 7 , with 7 ~ a s . Thus, twist-2 
DGLAP evolution is unable to describe the DDIS data 
below Q 2 = 5 GeV 2 at low £. We shall show that the Q 2 
and £ dependencies of the data arc, however, consistent 
with large contributions of HT effects. 

3. Twist decomposition of the saturation model. 

In the large energy limit the j*p scattering may be de- 
scribed in a colour dipole model [H [7] , where the 7* scat- 
tering is factorised into an amplitude of 7* partonic fluc- 
tuations and scattering of these states by multiple gluon 
exchanges. For DDIS one needs to consider the photon 
fluctuations into a colour singlet qq pair (a colour dipole) 
and into qq-ghion triple, that spans two colour dipoles in 
the large N c limit [1] . This gives the ^-integrated 7* cross- 
sections, dal"£/dM x = daf T /dM x +da q L %/dM x . As- 
suming an exponential ^-dependence of diffractive cross- 
sections, one finds for the qq component of 7*: 
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where bo is a diffractive slope and the sum over spins 
runs over massless (anti)quark helicities (h)h in the case 
of longitudinal photons whereas for transverse photons 



there is an additional average over initial photon polar- 
isations e, and / runs over the three light flavours. We 
use the photon wave functions ip^ {Q, z , r) m the form 
given in [5], and the GBW parametrisation [3] for the 
dipole-proton cross section <j(r) = a N(r) with N(r) — 
1 — exp(— r 2 /4i? 2 ). The saturation radius in DDIS de- 
pends on f , R(£) = (£/x a ) x/2 GeV- 1 and a = 23.03 mb, 
A = 0.288, x = 3.04-10~ 4 . 

We performed the twist decomposition of Eq. pi) 
through the Taylor expansion in inverse powers of QR. 

The contribution of the qqg component of 7* is calcu- 
lated at (3 = and in the soft gluon approximation (the 
longitudinal momentum carried by a gluon is much lower 
then carried by the qq pair). The correct /3-dependence is 
then restored using a method described in [6] , with kinc- 
matically accurate calculations of Ref. [5] . The soft gluon 
approximation is valid in the crucial region of M\ 3> Q 2 
or j3 <C 1, where the deviations from DGLAP are ob- 
served. With these approximations one obtains, in con- 
sistence with [5]: 
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where A = N c a 2 a s / 32Tr 3 b D and JVy = N(fj - fj), 
^01:^02,^12 = ^02 — ^01 denote the relative positions of 
quark and antiquark (01), quark and gluon (02) in the 
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to the dipole splitting kernel [TU] . The form of N 2 fol- 
lows from the Good- Walker picture of the diffractive dis- 
sociation of the photon [TT]. The factor \jM\ is a rem- 
nant of the phase space integration under the soft gluon 
assumption. The twist decomposition of ([3]) is performed 
using the Mellin transform in the r m variable: 
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where the expression for Hl,t(s) — J2f ^lt( s ) can 
be found in [5]. The Mellin transform of N 2 ^ „(roi) 
can be done in two steps. First one defines new in- 
tegrals N 2 q (s) = 1 1 — I 2 where, for a = 1,2 one has 
I a = (Q 2s M J dVoitrgi)- 1 / d 2 r 02 K 01l2 S a with St = 
(^02 + ^12 - N m N l2 ) 2 - {N m f and S 2 = 2N 01 {N Q2 + 
N12 — Nq 2 N\ 2 — Nqx). The integral 1% can be performed 
exactly, 

I 1 =tt (Q R) 2s 2 1+s (2 1+s - 1) T(s) 

s 1 
X [ff,- 3^2(1, 1,1-*; 2,2; -1)«], H s = J2t' ( 5 ) 

fe=l 
and for I 2 we use the large daughter dipole approximation 
ro2 > r 01 ,r 12 m r 02 and obtain, 
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FIG. 2. The LRG ZEUS data for &r w at low Q A compared 
to a DGLAP fit [3 and the DGLAP fit with included twist- 
4 and twist-4 and 6 corrections from the MMS saturation 
model. In yellow (gray) — the region of /3 where the correction 
due to qqgg may be neglected. 



These formulae agree within a few percent with exact 
numerical calculation of the diffractive cross-section |3| . 
The twist decomposition follows from Q as a sum over 
residues at the negative integer values of s. 

In the region of low £ the dominant contribution to 
the diffractive cross-section comes from qqg. However, 
at very low /?, an even more resolved photon fluctuation, 
with two emitted gluons, qqgg, may become relevant. 
Hence, at the lowest f3, one expects some underestima- 
tion of the dipole model predictions with qq and qqg 
states only. Therefore we use the data with /? > 0.035 in 
quantitative analysis. 

4. The evidence of higher twists observation. The 

twist decomposition of the saturation model provides an 
efficient tool to study the HT effects in data. We focus 
on the region of Q 2 < 6 GeV 2 , where the accuracy of 
DGLAP description breaks down. We compare the data 
to four different theory descriptions: (a) ZEUS-SJ NLO 
DGLAP fit [2] to all data with Q 2 > 5 GeV 2 , (b) the 
complete saturation model in the formulation of Rcfs. [6, 
[TT] (MMS-Sat), (c) the GBW saturation model @], and 
(d) the saturation models in which we retain only twist-2, 
twist-2 and 4 or twist-2, 4 and 6 contributions. In Fig. 2 
we compare selected results with data: the extrapolated 
DGLAP results, DGLAP plus MSS-Sat twist-4, DGLAP 
plus twist-4 and twist-6 MMS-Sat, and the other results 
are only described in the text. 

The saturation model results are obtained including 
qq and qqg components, using the original GBW param- 
eters A and a®, and three massless quark flavours. In 
our approach we modified the GBW parameter xq to 
£o = 2xo in order to account for the difference between 
Bjorken x and £, the variables used as V in GBW dipole 
cross-section in DIS and DDIS respectively. The variable 
£ corresponds to the actual momentum fraction flowing 



from the proton to the photon, and Bjorken x in DIS is 
only the threshold value of x in hadroproduction. The 
actual x of the gluon is larger, since typically, a hadronic 
mass M 2 ~ Q 2 is produced in DIS. We chose a s = 0.4 
that provides a good description of data. 

The conclusions from the analysis and from the figure 
are the following: (i) at twist-2 the DGLAP fit extrapo- 
lation, and the twist-2 components of the MMS-Sat and 
GBW agree with one another, but all fail to describe the 
data below Q 2 = 5 GeV 2 and at low £; the consistence of 
all approaches at twist-2 verifies reliability of saturation 
models at twist-2 (ii) the complete saturation models, 
MMS-Sat and GBW give similar predictions, better than 
the DGLAP description of the low Q 2 data, but still the 
predictions are much below the data; (iii) the combina- 
tion of twist-2 DGLAP and twist 4 part of MMS-Sat gives 
similar results as the complete MMS-Sat — the data are 
not well described; (iv) a combination of the DGLAP fit 
and twist-4 and twist-6 components of MMS-Sat gives a 
good description of the data at low Q 2 ; (v) in GBW one 
finds only twist-2 and twist-4 contributions, both close 
to those of MMS-Sat, and GBW is unable to describe 
the data with any truncation of the twist series. We also 
found that the large HT effects in the saturation models 
appear only in the qqg component. 

The quality of the improved description of the data 
with our best model, NLO DGLAP + twist-4 and twist- 
6 components of MMS-Sat is quantified with x 2 /d.o.f. 
in Fig 1, the dashed curve. The ZEUS LRG data [3] 
were fitted in this model for Q 2 > Q 2 ^^ In the fits only 
the DGLAP parameters were adjusted, and the twist-4 
and 6 components were taken from the MMS-Sat model, 
with fixed parameters a s = 0.4, £o = 2xo- Inclusion 
of HT terms greatly improves the fit quality in the low 
Q 2 region, — the maximal value of x 2 /d.o.f. ~ 1.5 at 
Q 2 nin — 2 GeV 2 is much lower than x 2 /d.o.f . ~ 3 of the 
DGLAP fit. The improved \ 2 is n °t Y e t perfect, but it 
should be expected, as the MMS-Sat model provides only 
a first and crude estimate of HT effects in DDIS. 

The final conclusions follow: (1) the discrepancy 
between the twist-2 DGLAP description and the data 
has strong Q 2 and ^-dependence, consistent with large 
contributions of twist-4 and twist-6; (2) the data require 
a truncation of the twist series in the saturation model. 

5. Discussion. We used an arbitrary truncation 
of higher terms in the twist expansion of the satura- 
tion models in order to determine a QCD motivated 
parametrisation of the £ and Q 2 dependence of the 
twist-2 DGLAP breaking effects in DDIS data. A trunca- 
tion of this kind, however, may be motivated in QCD. Let 
us recall that the GBW model assumes an eikonal scat- 
tering of the colour dipole off proton. This means, that 
one couples to the dipole an arbitrary number t-channcl 
gluons, without significant suppression. This is, however, 
in conflict with results of the BFKL framework (see e.g. 



Ref. [H]). In BFKL, at the leading logarithmic approx- 
imation, only one reggeized gluon may couple to a fun- 
damental colour line. Since DGLAP and BFKL approx- 
imations have the same double logarithmic (lnx InQ 2 ) 
limit, one concludes that also in DGLAP couplings of 
more than two gluons to a colour dipole is much weaker 
than in the eikonal picture. This, in turn, has implica- 
tions for the HT representation in the saturation model, 
because a n-gluon exchange provides the leading twist-n 
contribution in DIS cross-sections. Thus, we conclude 
that the eikonal saturation model should overestimate 
all the HT terms in single colour dipole scattering. 

In DDIS, however, also the qqg component is impor- 
tant, equivalent to two colour dipoles in the planar limit. 
Thus, one may couple up to four gluons (two per dipole) 
to the scattering state, without the BFKL constraint of 
one-per-line reggeized gluon coupling. Such a coupling, 
at the amplitude level, means that one may expect unsup- 
pressed contributions up to twist-8 in the cross-section, 
and a suppression beyond twist-8. This qualitative argu- 
ment provides possible motivation for truncations of the 
twist series in the saturation model, but in fact, in order 
to obtain reliable results, multi-gluon couplings to a qqg 
system need to be computed in QCD. 

Finally, let us point out that so far the HT terms 
in the saturation model were too weak to be seen in 
data, so a truncation of HT terms in saturation models 
should not lead to deterioration of other data description. 

Acknowledgements. The work is supported by 
the Polish National Science Centre grant no. DEC- 
2011/01/B/ST2/03643. 



* leszek.niotyka@u_i.edu.pl 

t mariusz.sadzikowski@uj.edu.pl 

* wojtek.slominski@uj.edu.pl 

[1] J. Bartels, K. J. Golec-Biernat and K. Peters, Eur. Phys. 

J. C 17, 121 (2000); J. Bartels, K. Golec-Biernat and 

L. Motyka, Phys. Rev. D 81, 054017 (2010). 
[2] J. C. Collins, Phys. Rev. D 57, 3051 (1998). 
[3] S. Chekanov et al. [ZEUS Collaboration], Nucl. Phys. B 

831, 1 (2010); Nucl. Phys. B 816, 1 (2009) 
[4] K. J. Golec-Biernat and M. Wiisthoff, Phys. Rev. D 59, 

014017 (1998), Phys. Rev. D 60, 114023 (1999). 
[5] H. Kowalski, L. Motyka and G. Watt, Phys. Rev. D 74, 

074016 (2006). 
[6] C. Marquet, Phys. Rev. D 76, 094017 (2007). 
[7] N. N. Nikolaev and B. G. Zakharov, Z. Phys. C 49, 607 

(1991). 
[8] M. Wiisthoff, Phys. Rev. D 56, 4311 (1997). 
[9] Y. V. Kovchegov and E. Levin, Nucl. Phys. B 577, 221 

(2000). 
[10] A. H. Mueller, Nucl. Phys. B 415, 373 (1994). 
[11] S. Munier and A. Shoshi, Phys. Rev. D 69, 074022 

(2004). 
[12] L. N. Lipatov, Phys. Rept. 286, 131 (1997). 



